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The appearance of the Fermi arcs or gapless regions at the nodes of the Fermi surface just above 
the critical temperature is described through self-consistent calculations in an electronic disordered 
medium. We develop a model for cuprate superconductors based on an array of Josephson junctions 
formed by grains of inhomogeneous electronic density derived from a pheise separation transition. 
This approach provides physical insights to the most important properties of these materials like 
the pseudogap phase as forming by the onset of local (intragrain) superconducting amphtudes and 
the zero resistivity critical temperature Tc due to phase coherence activated by Josephson coupling. 
The formation of the Fermi arcs and the dichotomy in k-space follows from the direction dependence 
of the junctions tunneling current on the d-wave symmetry on the Cu02 planes. We show that this 
semi-phenomenological approach reproduces also the main future of the cuprates phase diagram. 

PACS numbers: 74.20.Mn, 74.25.Dw, 74.62.En, 74.81.-g 



A convincing explanation to the origin of the pseudo- 
gap and the dome- like shape of the superconducting crit- 
ical temperature T^ip) on the doping level p of the copper 
oxide superconductors has become a major long standing 
challenge in condensed matter physics. The lack of a well 
accepted theory is likely to be due to the nanoscale com- 
plexity and the intrinsically inhomogeneous electronic 
structures [l|, that are not easy to be mathematically 
characterized. Recently, new techniques were developed 
to control the level of phase separation of dopant oxygen 
interstitials in La2Cu04+y and established a direct cor- 
respondence between the degree of order and the values 
of Ted, Q paving the way to new theoretical treatments. 

Electronic inhomogeneities in cuprates were verified by 
several different experiments: neutron diffraction [5(-[7|, 
muon spin relaxation (/iS'i?)[l, Q, nuclear quadrupole 
resonance (NQR) and nuclear magnetic resonance 
(NMR)[l3, have detected some form of disordered 
local electronic densities. The nanometer spatial varia- 
tions of the electronic gap amplitude A(r) measured by 
atomically resolved spectros copy such as scanning tun- 
neling microscopy fSTM) p3 - [la | is possibly related with 
the charge inhomogeneities. Angle resolved photon emis- 
sion spectroscopy (ARPES) found a large anisotropy in 
k-space: a larger gap at the leading edge of the Fermi 
surface or antinodal (along the Cu — O bonds) direction 
((±7r,0) and (0, ±7r))[l^ that remains well above and 
a small d-wave like gap that vanishes at the nodal direc- 
tion. 

More recently, ARPES experiments measured nodal 
gapless Fermi arcs starting at Tc and increasing with 
the temperature while the antinodal gaps remained al- 
most constants Another group revealed a metal- 
lic dispersion along the Fermi arcs and a Bogoliubov- 
like dispersion toward the antinodal regions either above 
and below Tcll^ demonstrating the presence of the su- 
perconducting amplitudes above Tc. Electronic Raman 



scattering experiments [l9[ showed that the nodal gap is 
connected with Tc and the antinodal with the pseudo- 
gap. Intrinsic tunneling spectroscopy also measured a 
gap that closes at Tc and a pseudogap that virtually 
does not change with the temperature and remains above 
r^(20j. in agreement with two energy scales in cuprate 
superconductors [TtI [l9|. 

To describe these complex phenomena some theories 
producing phase separation have been suggested, mainly 
based on doped Mott- Hubbard insulators |2l| - [25j . While 
they describe some of the observed features, they fail to 
predict all the details related with real space inhomo- 
geneities and the k-space dichotomy. While these theories 
describe some of the observed features of cuprates, they 
fail, for instance, to predict all the details related with 
real space inhomogeneities and the k-space dichotomy. 
To deal with this problems some theories considered 
the presence of impurities or inhomogeneities in the lo- 
cal electronic distribution as a starting point to derive 
the superconducting properties [2^ [13]. Some others ap- 
proaches used the influence of a mesoscopic phase sepa- 
ration in the appearance of superconductivity to derive 
the general properties of cuprates j28l[29|. 

In this letter we show that it is possible to find an uni- 
fied explanation to all these experiments starting with an 
electronic phase separation EPS forming small granular 
regions where isolated superconducting amplitudes may 
be formed. The origin of the EPS may be due to the 
lower free energy of the low density anti-ferromagnetic 
domains [10] ■ These local d-wave superconducting ampli- 
tudes in these islands form an array of multiple Josephson 
i unctions [31I [s^ . In a typical d-wave superconductors 
junction the direction dependence of the tunnel matrix 
elements that describe the barrier is relevant [s^- In the 
case of cuprates, all the domains have d-wave pair wave 
functions with the same direction with respect to the a 
and b crystal axis. Phase fluctuation above Tc generates 
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the gapless region that increases with temperatures at 
the nodal direction where the d-wave superconducting 
gap is smaller and accounts for the k-space anisotropy . 

To describe mathematically the electronic phase sepa- 
ration of high Tc oxides we use the time dependent Cahn- 
Hilliard (CH) equation [3^. In this approach it is possi- 
ble to follow the formation of domains with distinct local 
densities below the phase separation transition tempera- 
ture Tps- Tps ~ Tps{p) is likely to be close the upper 
pseudogap temperature since Tps{p) may be the cause 
of some observed high temperature anomalies [ssf. The 
time t in the CH equation is related with the temperature 
T; larger times correspond to lower temperatures below 
Tps(p) and consequently, larger disorder j29l Isol IbgI. . 
In the disordered phase the appropriate order parameter 
is the difference between the local and the average charge 
density p, u{p,i,T) = {p{i,T) —p)/p. u{i,T) ~ corre- 
sponds to the homogeneous case above Tps{p) and the 
Ginzburg-Landau (GL) free energy functional is given by 
the usual power expansion, 




FIG. 1. (color online) Typical solution of the free energy 
potential Vgl{u, T) on a 50 x 50 lattice. At low temperatures 
the charges are attracted to the grains, similar to a granular 
superconductor. On the top left is shown the d-wave super- 
conducting amplitude in these bound regions. 
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+ Vgl{u,T). 
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Where the potential Vgl{u,T) = -A'^{T)u'^/2 + 
+ A^{T)/B = a{Tps{p) - T) for T < Tps, a 
is a constant, e gives the size of the barrier between the 
low and high density phases [2^, . The CH equation 
can be written [ssj in the form of a continuity equation 
of the local density of free energy /, dtu = —V.J, with 
the current J = MV{6f/6u), where M is the mobility 
or the charge transport coefScient that sets the phase 
separation time scale. Therefore, 



du 
'dt 



-My^{s^\7^u- A\T)u + B^u-^). (2) 



In Fig.(IT]), we show a typical VGLiu-,T) simulation 
where the barriers between different grains are clearly 
visible. We have argued in previous works [Sll [32| that 
the particles may become confined in these domains in- 
creasing the probability of Cooper pair formation. In our 
calculations we take into account that these free energy 
barriers are essentially constant at low temperatures and 
varies as (Tps — T)^'^ near the transition as shown by 

chUI. 

Then, after performing the CH simulations, the lo- 
cal disorder density p{i,T) is used as the initial input 
and it is maintained fixed throughout the self-consistent 
Bogoliubov-deGennes (BdG) method. This method 
yields local varying d-wave amplitudes Ac{i,p,T) as 

demonstrated previously [13-111: 113 ■ 

As mentioned above, Vgl{p,T) favors the charge to 
be confined into the domains acting as a catalyst to 
Cooper pair formation. We have showed that the two- 
body attractive potential of the BdG equation is scaled 
by the height of the free energy barriers shown in Fig.([3|). 
In Ref ■ ( |32| ) we called it the grain boundary poten- 
tial Vgb- Here we generalize it to any doping value by 



the expression, Vgb{p) = -14.2 + 45 x p (eV). With 
this two-body attractive potential a typical BdG solu- 
tion Ad{p = 0.15, z,T) is shown in Fig. ([2]) for five dif- 
ferent grains "i". The low temperature values between 
30-70meV are in the energy range of the LDOS gaps mea- 
sured by STM experiments on Bi2212[T2i^. 



70 
60 
50 

I 40 
S 30 
20 
10 





site 1 




site 2 




site 3 ' 




■^---^ site 4 




site 5 ' 




Tc \ ' 



40 



60 



80 



100 
T(K) 



120 140 160 



FIG. 2. The temperature evolution of Ad{i) at five different 
locations of ap = 0.15 compound. The mean-field BdG calcu- 
lations all vanishes at the same temperature T*(0.15) « 165K 
(above rc(0.15) — 92K). The low temperature values between 
30-70meV are in the energy range of the Bi2212 LDOS gaps 
measured by STM. 

To obtain the experimental values of Tc{p), we use 
an approach similar _to a system composed of gran- 

Our main proposal is that 



ular superconductors [41 



the superconduc ting transition occurs in two steps as T 
decreases [Sli [3^ l37j : first by intra-grain superconductiv- 
ity and than by Josephson coupling with phase locking 
at low temperatures. 
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These two completely different calculations, yield- 
ing two different energy scales, are motivated by the 
two ener gy scales found in most cuprates mentioned 
above [T6lTl9|. and also in the two different regimes of the 
fluctuation magnetoconductivity in a Y Ba2Cu30'j sin- 
gle crystalHij. These systematic measurements detected 
an effectively two-dimensional (2D) regime far above Tc, 
that can be interpreted in our approach by the isolated 
superconducting grains in the CuO without phase coher- 
ence. Decreasing the temperature towards a crossover 
to a three-dimensional (3D) Gaussian regime sets in as 
the different superconducting regions or grain become 
connected and increase coherence. Very close to a 
critical regime characteristic of a 3D XY university class 
is measured which is consistent with a granular supercon- 
ductors in which each superconducting grain develops its 
own phase $ that may oscillate but become locked at Tc, 
yielding long range order. 




FIG. 3. Schematic figure representing the array of Josephson 
junctions on a CuO'z plane. The local dependent amplitudes 
A(p, i, T) are represented by the orbitals and the phase angles 
cj> oscillates around the crystal a-axis direction due to the 
symmetry of the d-wave on the CuOz plane. 

Consequently, the system is regarded as an array of 
Josephson junction as shown schematically in Fig. ([3]). 
The phenomenology of Josephson tunnel junctions be- 
tween d-wave superconductors have been developed by 
Bruder et aljs^. They calculate the tunnel matrix ele- 
ments in second-order perturbation theory for two su- 
perconductors (1 and 2) with superconducting ampli- 
tude Ad,i/2(i,r,$) = A(i,T)cos[2($ - 01/2)]. We use 
their results to an array of junction as schematically 
shown in Fig. ([3]), where the local amplitudes Ad{i,T) 
are represented by the size of orbitals and the phase 
angles (j) have been drawn around x-direction. $ is 
the polar angle of the first Brillouin zone. Due to the 
thermal energy the gap must be averaged over (j), then 
< A(p,T,$) (/'^*A3''(p,T)cos[2($-</>)]d</.)/A</.. 

And A7'(P, r, ) = Ad(p, i, T)/N is the spatial aver- 
age superconducting amplitude. At T < Tc the values 
of 01/2 are essentially locked. At T > Tc the phase 
angles decouple and oscillate within a maximum am- 
plitude that increases with T as shown in Fig. ([3]). In 
this case < A™(p, T)cos[2($ — (j))] > vanishes when 



2($ -</))« ±7r/4 or it is less than KbT. The resuhs 
are shown in Fig.Q and compare well with the $ de- 
pendent ARPES measurements of Lee at al[l3|. The ap- 
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FIG. 4. The absolute value of the spatial and (j> average, 
< A2''(T, $) >, as function of the k-space polar angle <I>. 
The developing of the Fermi arcs around the nodal direction 
above Tc is a consequence of the lack of phase coherence and, 
as measured by ARPES [13, El , < A7(r,<I>) >< KgT. 



pearance of a gapless region around the nodal direction 
($ = ±7r/4) above Tc while the antinodal ($ = ±ri7r/2, 
n=l,2,..) gap remains almost constant is concordant with 
the data of Kanigel et alp^]. They found that the elec- 
tronic dispersion is either metallic along the Fermi arcs or 
exhibits Bogoliubov-like behavior near the nodes at tem- 
peratures above Tc. This behavior can be understood 
with the results for Ad{i,p,T) x T shown in Fig. ([2]) and 
A7(p,T) X $ of Fig. (31). 
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To calculate the Josephson coupling energy that is con- 
nected with the onset of phase coherence, we recall the 
work of Bruder et aljs^l- They found that the tunneling 
current behaves in a similar fashion of s-wave supercon- 
ductors junction and the leading behavior is determined 
by tunneling from a gap node in one side of a junction into 
the effective gap in the other side. Consequently, as a first 
approximation to the Josephson coupling energy Ej, we 
adapt the theory of s-wave granular superconductors [isj 



to the A7(p,r) in the grains. 



Ejip,T) 



^ 2KrT ' 



4e2i?„(p) 



(3) 



Where i?„(p) is the normal resistance of the compound. 
We use here values proportional to the planar resistiv- 
ity Pab measured on the La2-pSrpCu02 series[4J] and 



rpCu(J2 series 
the estimated values of Rn{p)e^ /h are plotted in Fig.©. 
Since the onset of phase coherence occurs when the ther- 
mal energy equals the Josephson coupling, KbT{p) = 
Ej{j),T), we obtain one of our most important result; 
the derivation of the dome-like shape dependence of 
Tcij}) = Ej{p,T,)/Kb, as plotted in Fig.®. 
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FIG. 5. The temperatures which Ej{p,T) = KbT{p) yield 
Tc{p) = Ej(p,Tc) / Kb, for comparison we plot also the for- 
mula rc(p) = 95(1 - 82.6(p - 0.16)^). The experimental re- 
sistivity Rn{p)e^/h used in Ej(p,T) (Eq|3} is also plotted 
following Takagi et aljZJ- It is also shown the average low 
temperature gap A^" (p, T = 0) that is proportional to the 
pseudogap temperature T*{p). 



To show how Tc{p) is calculated we plot also in Fig.® 
the average low temperature gap A'^"{p,T — 0) 
A°i^{p,Tc) that enters into Eq.Q. As it is shown in 
Fig.®, the results agree with empirical formula Tc{p) = 
95(1 — 82.6(p — 0.16)^). The dome-like shape is a con- 
sequence of the values of AJ^" {p, T) that increases as p 
decreases while l/i?„(p) have opposite behavior, it is es- 
sentially zero near p = 0.05 and increases with p (see 
Rn{p)e^/h plotted in Fig.dl])). 



In conclusion, the CH disordered free energy in 
cupratcs forms small electronic domains that induces 
local Cooper pairs and also an array of Josephson 
junctions, as in a granular superconductor. This ap- 
proach reproduces important experimental results and 
provides novel interpretations: The two e nerg y scales 
measured by electronic Raman spectroscopy [19[, tunnel- 
ing spectroscopy [13 and ARPES[l^[l3 are; one associ- 
ated with the Josephson coupling energy Ej{p) that fol- 
lows Tc{p) and the other related with the d-wavc mean 
superconducting gap A|^" (p, T) that increases as the dop- 
ing decreases. 



Concerning the k-space anisotropy, the appearance of 
the Fermi arcs near the nodal direction is due to the lack 
of phase rigidity among superconducting domains com- 
bined with the small d-wave amplitudes around the nodal 
directions. In this way the existence of Bogoliubov-like 
dispersions above Tc, close to the antinodal directions, 
have the same cause, namely the anisotropic d-wave ge- 
ometry of the intragrain superconductivity that is larger 
around these directions and remains just at T*{p), i.e., 
above Tc{p). Consequently the pseudogap phase is com- 
posed of these localized superconducting domains with- 
out phase coherence and with a nodal metal structure. 
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